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Abstract 

A new approximation scheme for non-perturbative calculations in a quan¬ 
tum field theory is proposed. The scheme is based on investigation of solutions 
of the Schwinger equation with its singular character taken into account. As a 
necessary supplementary boundary condition the Green functions’ connected 
structures correspondence principle is used. Besides the usual perturbation 
theory expansion which is always available as a particular solution of our 
scheme some non-perturbative solutions of an equation for the propagator 
are found in the model of a self-interacting scalar field. 

Introduction 

As is well known there exists a wide class of physical phenomena (called non 
perturbative effects) which cannot be described by a finite number of terms of the 
perturbation theory series. It is usually believed that the full or partial summation 
of the perturbative series solves this problem. Such an approach supposes implicitly 
that the sum of the perturbative series contains exhaustive information about a given 
quantum field model. 

Meanwhile even simple mechanical models provide an example for the perturba¬ 
tive series to be unable to describe a physical situation properly. 

Thus, it is well known from fluid mechanics that one cannot use an ideal liquid 
as the leading approximation to a viscous one (even if the viscosity is very small) 
when considering a process of flowing near the boundary of an immersed body. It is 
a consequence of the fact that the viscous liquid is a so called singularly perturbed 
system 11 compared to the ideal liquid. A solution of equations of motion for the 
ideal liquid as well as the perturbation theory based on this solution cannot obey 
the boundary conditions for the viscous liquid in principle. 

A characteristic feature of a singularly perturbed system is existence of an es¬ 
sential singularity with respect to an expanding parameter in solutions. If such a 
system is described by a differential equation then the highest derivative coefficient 
depends on an expanding parameter as follows: when this parameter is equal to zero 
the order of the equation decreases and the number of necessary boundary conditions 
reduces. This means that the perturbative series is a solution of the initial problem 
only for some special boundary conditions. If a physical situation goes beyond their 
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framework the result of the perturbation theory summation fails to describe such a 
situation in principle. 

Let us consider an elementary example. Take the equation for a system in the 
form: 

A x = t — x (1) 

with the boundary condition x(0) — X. At A = 0 the order of the equation lowers 
so ([[]) is a singularly perturbed equation. In this case the perturbative series on A 
breaks at the second term and the result of its summation is 

Xpert (t) = ~ ( 2 ) 

Though x per t is an exact solution of (P with the special (perturbative) boundary 
condition x per t( 0) = —A, it cannot serve as a solution of the problem (P when 
X —A. The exact solution of ([I]) has the form 

x(t) = t - A + (X + A) exp(— y) (3) 

A 

and shows that in solving a problem with non-perturbative boundary conditions 
by iterations one should introduce a so called boundary series besides the usual 
perturbative one. In our simple example (|l|) this boundary series consists of the two 
terms containing an essential singularity on A. 

The boundary series contribution cannot be taken into account by perturbation 
theory and in general this contribution dominates near the point t — 0 that is 
outside the region of the Tikhonov’s theorem applicability 0. This domination can 
be easily seen from our example if we calculate the derivatives of x(t) at t = 0. We 
get i(0) = —\X, x(0) = AX + A and so on. These derivatives are of obviously non- 
perturbative character (except for the case of the perturbative boundary condition 
X = —A). 

Now let us turn to the case of quantum field theory. An essential singularity with 
respect to a coupling constant is the wide known attribute of any field theory model 
with interaction. Moreover, the Schwinger equation for the generating functional 
contains a coupling constant as a coefficient at the highest functional derivative so 
from the point of view of the differential equations theory the Schwinger equation is 
a typical example of a singularly perturbed equation. 

It is also worth mentioning that it is Green functions (vacuum expectation values) 
which are of physical interest. As these functions are the derivatives of the gener¬ 
ating functional at the source switched off (hence, in the obviously non-Tikhonov 
region) so the boundary series contribution to Green functions will dominate ex¬ 
cept for the perturbative boundary conditions. In this work we make an attempt 
to treat the Schwinger equation as singularly perturbed one and to construct a new 
approximation scheme of solving of this equation. 

The proposed scheme approximates the non-perturbative Green functions already 
at first step without summing of expansions and at the same time contains the stan¬ 
dard perturbation theory as a particular case. The scope of admissible boundary 
conditions becomes wider and they do not limited by the framework of the per¬ 
turbation theory any more. It opens new possibilities for describing of essentially 
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non-perturbative effects (such as mass spectrum reconstruction) already at the lead¬ 
ing orders of approximation. We hope that the proposed approximation scheme (or 
those similar to it) can prove to be useful when studying the problems connected 
with the spontaneous symmetry breaking and also in problems where one should 
take into account from the very beginning the nontrivial properties of the physical 
vacuum as media. 

Our work is organized as follows. In part 1 necessary definitions are introduced 
and general principles of the approximation scheme construction are stated. The 
primary version of the scheme || based on a general solution of the Schwinger equa¬ 
tion is considerably simplified and in the form proposed here this scheme can be 
generalized to any quantum held model (in present work we deal with the simplest 
case of a self-interacting scalar held only). In part 2 the equation for the propagator 
at the leading approximation is solved. Brief discussion of the results contains in 
Conclusion. 


1. Schwinger Equation and Approximation Scheme 

We will consider the theory of complex scalar held (j>(x) in d -dimensional Eu¬ 
clidean space x € E d with nonlocal self-interaction: 

Sm = J dx dy {4>*(f>) (x) ( 0 * 0 ) (?/). (4) 

The limit A (xy) —> \S(x — y) corresponds to the local theory with the quadric 
interaction. 

To define the Green functions we introduce a bilocal source rj(xy)f _ ]. Then the 
n-particle ( 2 n-point) functions (vacuum expectation values) are the derivatives of 
the generating functional 
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The hrst derivative is the propagator of the particle: 


A (xy) = G i 
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( 6 ) 
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the second one is the two-particle function: 
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1 The use of the bilocal source is made only for the sake of simplicity and compactness of calcu¬ 

lations and is not a principal point in proposed formulation of the approximation scheme. 
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The Schwinger equation for the generating functional is a consequence of the held 
equations and the quantization conditions. For the model involved it has the form: 

S 2 G 


dx' A (xx') 


[m 


d [ 


S7](yx') dy^x'x) 

dG f / . dG . . 

- + dx n[xx ) ———— + d(x — y) G. 
) J dr]{yx') 


5rj{yx 


( 8 ) 


Here m 2 stands for the mass of the free held cj)(x) and d 2 = J2 n =i ~§x~ is the Laplace 
operator. 

The following notations will be subsequently used: 


8 (5 

(AB)(xy) = A(xy)B(xy), = 

(A * B) (xy) = J dx'A(xx') B(x'y). 

In these notations the Schwinger equation (IS]) reads as follows: 


dG 
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dG 


A -r— * —— (m 2 — d + r/) * —-h G 
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( 9 ) 


( 10 ) 


Successive differentiation of m gives (on switching oh the source) the inhnite 
system of the Dyson equations for the Green functions G n . 

The Schwinger equation fllO|) belongs to the class of singularly perturbed equa¬ 
tions. Indeed, when A = 0 the order of the equation reduces: the generating func¬ 
tional CA 0 ) at the leading order of the perturbation theory on A (that is the free held 
generating functional) is a solution of the hrst order equation: 

6G<® 


(m 2 — d 2 + r )) * 


drj 


+ G (0) = 0. 


( 11 ) 


The solution of (JTJl) is uniquely hxed by the only boundary condition: the normal¬ 
ization of the generating functional 


G [rj = 0] = 1. 

This solution is: 

G(°) — d e t | (m 2 — d 2 + 77 ) _1 * (m 2 — d 2 ) 
The perturbative solution for the Schwinger equation 

OO 

G per t = £ G( "’. 


( 12 ) 

(13) 

(14) 


n=0 


where G (n) = 0(A n ) is constructed by the iterations. The equation for G^ n ' is again 
of the hrst order with respect to the functional derivatives 
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and therefore to fix at any order of the perturbation theory on A one needs 
the only condition G^[ 0 ] = 0 which is a simple consequence of the normalization 
condition (|T 2 |). 


Meanwhile at A 7 ^ 0 (even if small) the Schwinger equation (|T 0 |) is that of the 
second order with respect to the functional derivatives and for its solution to be fixed 
uniquely one needs a supplementary boundary condition. The problem of bound¬ 
ary conditions is understood to be the key one for constructing of non-perturbative 
solutions of the Schwinger-Dyson equations IB SB’ 

Being the solution of the Schwinger equation fllOD , the perturbative series (|l4l) 
has the doubtless advantage that there is no need for fixing a supplementary bound¬ 
ary condition. But this advantage can turn to be a serious lack: physical phenomena 
described in principle by a given held model can escape out of the consideration. 
The perturbation theory even being summed up can fail to obey the physical bound¬ 
ary conditions just as the ideal liquid is unable to provide the physical boundary 
conditions for the viscous liquid in the boundary layer. At the same time enlarg¬ 
ing of the class of admissible boundary conditions can provide the description of 
non-perturbative phenomena after making a finite number of iterations and without 
summing of the corresponding expansions. 

Let us turn to the construction of the approximation scheme. Since the per¬ 
turbation theory is the only universal tool for calculating of the Green functions in 
quantum held theory we will use it as a base that is will construct our scheme in 
such a way that the perturbation theory expansion would be contained in it as a 
particular case. 

Introduce a perturbative approximant: 


N 


G", = £ G (n) 


(16) 


n =0 


and define a non-perturbative approximant G N to be found as a solution of the 
equation: 
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(17) 


When N —> 00 the perturbative approximant G^ ert tends to the perturbative solution 
of the Schwinger equation G pert and therefore G N tends to some exact solution of 
the same equation with, in general, arbitrary boundary conditions!- That is why it 
would be more correct to speak about a set of approximants {G^} every element of 
which obeying def ini te boundary conditions. Under N —► 00 this set of approximants 
turns into the full set of the all solutions of the Schwinger equation. G^ ert G {G^} 
as the trivial solution of equation (|T7|) . 

2 We do not discuss the problem of convergency of expansions. 
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To construct a solution of equation ([L7|) one needs to formulate some criteria 
restricting the framework of the admissible supplementary boundary conditions so 
as in general case they are arbitrary. After such a restriction we obtain a certain 
subset {G N } of the all possible solutions of equation (|l7|): {G N } C {G*^}. We will 
choose the boundary conditions in such a way (see below) that the set {G N } would 
include G pert . Our scheme will be nontrivial if it will contain other solutions too. In 
such a case one comes again to the question of the choice of the proper solution for 
the given model among the elements of {G N }- This choice can be made on taking 
into account some additional physical requirements. One of them is the principle of 
minimal energy value of the ground state. 

So, summing all mentioned above, our goal is to reduce the full set of the all possi¬ 
ble solutions {G^} to some subset {G N } by imposing additional boundary conditions 
on the solutions of equation (|T7|) . 

In formulating such conditions we will base ourselves on those properties of the 
perturbative approximant which can be generalized to the non-perturbative case. In 
so doing we extremely restrict the admissible boundary conditions but the pertur¬ 
bative solution is still contained in {G N }- 

One of the properties mentioned above is the connected structure of the pertur¬ 
bative approximant. General connected structure of Green functions can be found 
beyond the framework of the perturbation theory with the help of the theorem on 
the connectivity of the logarithm^. For example, the connected structure of the 
two-particle function is: 

= A(x lVl ) A (x 2 y 2 ) + A(x lV2 ) A (x 2 V i) +G ? n (*y l ylj • ( 18 > 

G 2 on being the connected part of the two-particle function and A is the full propa¬ 
gator. For the three-particle function this structure is of the form 

G 3 = Sym {6 A A A +9 A G c 2 m } + G c ° n (19) 

where the notation Sym stands for the Bose-synnnetrization. Formulae G-0) 
(and those similar to them for higher functions) are valid not only at any order of 
the perturbation theory but also beyond its frames. 

At N- th order of the perturbation theory all the connected parts of n-particle 
functions are equal to zero when n > N + 1. For example, at the leading order only 
the propagator is not equal to zero, at the order 0(A) (N = 1) the propagator and 
the connected part of the two-particle function are non-vanishing, and so on. 

To say it another, dynamics at IV-th approximation is fully defined by N + 1 
lowest Green functions. In this connection at IV-th order one can confine oneself to 
considering the subsystem of IV + 1 equations from the whole infinite Dyson system. 
So as at N —> 00 the N + 1 equations turn into the full system then the truncated 
approximant (a solution of N + 1 equations) still tends to the solution G pert . 

3 It is worth mentioning that to find the connected structure one should introduce the simple 
sources. The derivatives of Z = InG with respect to the bilocal source are not in general the 
connected functions. 
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We lay these properties into the base of our scheme construction. Physically this 
means that we consider the non-perturbative solutions of the Schwinger equation of 
such a type that the highest Green functions weakly affect dynamics defined by the 
lowest functions that is, for example, one can neglect the three-particle forces when 
considering the leading approximation to the two-particle processes. 

Thus, to construct the A-th step of the approximation scheme one should proceed 
in the following way: 

1. At the A-th step we consider A + 1 equations of the Dyson type following 
from ([L7|) that is equation ([17) itself and its A derivatives at g = 0. 

2. Due to A + 1 equations contain (A + 2) n-particle functions G± ,..., G^ +2 this 
system needs a supplementary boundary condition^. To introduce such a condition 
we will demand the correspondence between the connected structures of the sought 
approximant and the perturbative one. This perturbative approximant consists of 
the set of (A + 2) n-particle functions calculated with the perturbation theory which 
will be denoted as 


g N n = (- 1 ) 


8 n G 


N 

pert 


5 rj n 


77=0 


For the perturbative approximant [g^+ 2 ) c ° n = 0 therefore we put for our system 


r N 

^ n +2 


con 


0 . 


( 20 ) 


The connected structures correspondence condition ( 0 ) allows one to express 
(A + 2)-particle function in terms of A + 1 lowest ones 


7 TV 


— fN nN n 

— b_/\r+2 


7V+1 


( 21 ) 


and by this to close the system of equations. 

At every step of the approximation scheme we get, generally speaking, some 
collection of the n-particle functions sets {Q± ,...,Q^ +2 }. The set of perturbative 
n-particle functions {g^, ..., 9 ^+ 2 ) is a ls° contained in the collection. A priori the 
region of the scheme applicability is that of the small A. When A —> 00 each of 
these sets tends to the full (infinite) number of Green functions corresponding to 
some exact solution of the Schwinger equation (0)- Each of these solutions has the 
true connected structure. 

Let us elaborate the two primary steps in more detail. For A = 0 the perturbative 
approximant G pert = is the generating functional of the free Green functions 
(|I3D. The system consists of one equation — this is equation (p~7|) at r/ = 0: 


dx' A (xx 1 ) Q 01 X X ' 1 " 2 • a2 


2 \ yx , I +(m -d ) A u (xy) - S(x - y) 


= j dx' A (x x') g ° | 


yx! 


( 22 ) 


Normalization condition © is supposed to be fulfilled at any step of consideration. 
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Here 


and 


A N = - 


sg 
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Sr] 


Tj =0 


92 ( f. 1 X T 2 ) = A c (xij/i) A c (x 2 j/ 2 ) + A c (xi2/ 2 ) A c (x 2 yi) 

\Vl x 2 


(23) 


(24) 


- the two-particle function of the free held where A c = (m 2 — d 2 ) is the free 
propagator. In accordance with the connected structures correspondence condition 
(| 20 |) [G^) c ° n = 0 that is we have for the two-particle function Q 


gl ('V X . 2 j = A°(xi?/i) A 0 (x 2 y 2 ) + A°(xi y 2 ) A 0 (x 2 yi). 
\ y\ 1J2 j 


(25) 


So we have the equation for the propagator A 0 at the leading approximation. 

At the next step (N = 1 ) the system contains two equations. As it follows from 

(HD 


(m 2 — d 2 + rf) * 
so the system can be written as: 
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+ Gp ert ~ ( ^ 
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J dx'X (xx 1 ) + (m 2 — d 2 ) A 1 (. xy) — S(x — y) — 


dx'X (xx') 


( rn \ ( nr* rtr 

Jb JL \ Ql Jb Jb 
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(26) 


f dx'X ( Xl x') g}i XxX2X '\ + (m 2 - d 2 xi ) g{ Xx X2 
J {ymx' J xlJ \yiy 2 

<5(xi - y 2 ) A 1 (x 2 y 1 ) = J dx'X (xjA) 


S(x 1 - yi) A 1 (x 2 y 2 ) 


[f Xi X 2 X 

9 Am y2 x' 


of Xi X 2 X 

9 A yi 92 x’ 


(27) 


At N = 1 that is at the order 0(A) of the perturbation theory the connected part 
of the three-particle function $3 is equal to zero. Hence the connected structures 
correspondence condition at this step reads: 


Ql 


= 0 


(28) 


and in accordance with (|I8|-[I9[) equations (|26D and (p7| ) form the system for A 1 and 

[GIF”. 

To the end of this section we point out another way of the approximation scheme 
construction based on reducing the order of the Schwinger equation. 

Let us recall that if any particular solution of a linear n-th order differential 
equation is known the equation can be reduced to that of (n — l)-th order. This 
procedure can be applied to the Schwinger equation O treated as a linear second 
order differential equation with respect to the functional derivatives [7]. 
















Define a functional W in accordance with the formula 


W(xy\rj) = R ( xy\rj ) G[rj\ 


5G 

5r](yx) 


where R is a solution of the following equation: 


( 29 ) 


i'A Rf t — +(Afi) * R= 1 + (m 2 -d 2 + ri)*R. 
or] 


(30) 


The equation for W follows from ( |T0|) and ( |30|) and turns out to be of the first 
order: 


(XW) * — +(\R) *W = (m 2 - d 2 + y) *W. 
or] 


(31) 


To use this method one should find at least one exact solution of equation @ 
or equivalently the particular solution G part of ([H]) (then R par t = ^ In G part ). Such 
a solution was found in [6] (see also 00) and is of the form: 


Rpart = J dx'dy'X l {yy')r](y'x') +m 2 J dy'X 1 {yy'). (32) 


On the other hand one can perform “approximate” lowering of the order by 
using the perturbative approximant R N which is defined as being an approximate 
perturbative solution of (^0|) up to 0(X N ) accuracy. 

On solving the equation for W N 


(A W N ) * +(XR N ) *W N = (m 2 - 8 2 + 77 ) * W N (33) 

or] 

we obtain the sequence of approximants {R N , W N } and from (|29|) can recover all 
the derivatives at y = 0 of the corresponding functional G N . For the trivial solution 
hF A = 0 we get G N = G^ ert — the perturbative approximant (|16j) and for any 
nontrivial W N ^ 0 we hnd some non-perturbative approximation. The problem of 
supplementary boundary conditions can be solved with the help of the connected 
structures correspondence condition. The final equation for the propagator at the 
leading approximation coincides with (j22|). 

But it should be mentioned that this way is more technically complicated because 
of the necessity to take into account nontrivial integrability and Bose-symmetry 
conditions for W. An advantage of this scheme is the possibility to use as a starting 
point the exact solution (|32|) . But in so doing the question about the boundary 
condition should be solved anew because it is clear that R per t = hirijv^oo R N 7 ^ R par t- 


2. The solution of the equation for the propagator 

At the leading approximation (N = 0) equation (^) with ( |53| ) (^) gives the 
equation for the propagator A 0 which will be denoted as A in this section. 

Note, that at this step of the approximation we cannot use the principle of mini¬ 
mal energy value of the ground state in order to fix the physical solution because for 
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so doing one should investigate (as minimum) equations of the next step (N = 1) in 
which the connected two-particle function is contained. 

Nevertheless already at the leading approximation it is possible to solve the ques¬ 
tion of nontriviality of our scheme that is to verify if the equation for the propagator 
admits solutions different from A c . If exist, these nontrivial solutions will be of the 
nonperturbative origin by the construction. 

The equation for A can be written in the form: 


(m 2 — <9 2 + E) * A = (m 2 — d 2 + E c ) * A c 
The operators E and E c are formally defined as: 

E(xy) = A (xy) A (xy) + 5(x — y) J dx\X(xx\) A (aqaq), 

E c (xy) = A (xy) A c (xy) + S(x - y ) J dx 1 X(xx 1 ) A c (x 1 x 1 ), 


(34) 


(35) 

(36) 


The structure of these quantities coincide with that of the mass operator at the 
order 0(A) of the perturbation theory. In the local limit A (xy) —a XS(x-y) they are 
the formal divergent quantities and the problem arises of equation ([3f|) renormalizing. 
Due to E is expressed in terms of the propagator to be sought it is sufficient to 
define E c . The usual way of handling such a quantity is to introduce a regularization 
E c regT, c and then to define a renormalized quantity by subtractions. However 
this technique is inconvenient and complicated in our case. 

More simple recipe is to apply “the renormalization without subtraction” method 
[HI- This method as well as that of “differential renormalization” 10, |I0|| is one of 


the realizations of Bogolubov’s idea to define products of distributions without using 
senseless divergent quantities. The general idea of the method is to define such a 
product as a distribution from the Schwartz space being a solution of some equation 
well defined in this space. 

In the perturbation theory this method is equivalent to the usual ones. The 
results thus obtained can differ only by finite renormalization as they do for different 
subtraction procedures. 

The quantity E c to be defined consist of the two terms. At the local limit these 
terms are equal to each other and E c can be formally written as 2A A c (O)h(x — y). 
Due to A c (x) has a singularity in the origin when d > 2 this formal expression needs 
to be correctly defined. 

Define this as a solution of an equation in the Schwartz space. The equation can 
be found as follows. Let us consider the quantity: 


(x - y ) 2 A c (xy). 


In d < 4 this is regular at x = y because the singularity of the propagator A c is 
“killed” by the factor (x — y) 2 . Therefore the expression (x — y) 2 A c (xy)X(xy) is a 
well defined distribution at the local limit A (xy) —> X5(x — y). We lay this quite a 
definite in sense of distributions limit into the base of E c definition. Namely we will 
take E c as a solution of the equation 


(x - y) 2 Z c (xy) 


2A ((x 


yf A c (xy)) 


5(x 


y)- 
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i) d — 2. In the two-dimensional space we have under x —> y 


A c (xy) = In 


[x 


y) 2 m 2 


+ ip(l) + O ((x — y) 2 In (x — y) 2 rri 2 'j . (37) 


If follows from that at local limit A (xy) —> \5(x — y) the mass operator £ c 
satisfy the equation: 

{x ~ y) 2 E c (xy) = 0. (38) 

We will define £ c as a 0(2)-invariant solution of (|38|) belonging to the Schwartz space 
S / (E 2 ). Such a solution is 

T, c (xy) = C 5(x - y) (39) 

where C is a constant. It is a finite mass renormalization as can be expected. From 
formula ( |37| ) one can also find that the difference of the two free propagators with 
masses rri \ and m 2 

I 777,2 

A c (xy\mi) - A c (xy\m 2 ) = — In + 0((x - yf ln[(x - yfm 2 }) (40) 

is non-singular at coinciding arguments. Therefore at d = 2 

A m 2 

£ c (au/|mi) - E c (xy\m 2 ) = — 8(x - y) In —§. (41) 

Z7T 777-^ 

That is a hnite renormalization cancels in the difference of the two mass operators. 
It is a reflection of the fact that the divergencies of the corresponding integrals in 
the momentum space cancel. 


ii) d — 3. This case is similar to the previous one. In the three dimensional space: 

A c ( xy ) =- 1 exp{—mJ(x — y) 2 } (42) 

4 n^{x-y ) 2 

and £ c is 0(3)-invariant solution of (^) belonging to the Schwartz space S'(E 3 ) and 
is given by the same formula (^). The two mass operators difference also does not 
contain the renormalization arbitrariness and reads as follows 

£ c (au/|mi) - E c (xy\m 2 ) = 5(x - y) (m 2 - mi). (43) 

Z7T 

Here m = Vm 2 . 

iii) d — 4. In this case the theory contains quadratic divergencies and result 
changes. The free propagator is more singular at coinciding arguments 

a c (xy) = 4?r2 , s 2 + QCMfo - v) 2 ™ 2 ])- ( 44 ) 
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As a consequence we get an inhomogeneous equation for £ c at the local limit 

{x-y) 2 E c (xy) = ^8(x-y). (45) 

Its 0(4)-invariant solution from the Schwartz space S^E^) is 

K{xy)= iL^ d 25{x ~ y)+c5{x ~ y) (46) 

that is £ c contains the finite term of the wave function renormalization besides the 
mass renormalization. This term reflects the presence of the quadratic divergencies. 
In contrast to the previous cases the difference of two £ c 

E c (xy\mi) - E c (xy\m 2 ) = C'5(x - y) (47) 

contains a finite arbitrariness. To say another, in the corresponding momentum 
integrals the quadratic divergency cancels but the logarithmic one still remains. 

Having defined £ c in ( |36| ) as described above let us turn to the solutions of 
equation (|34|). In what follows we will everywhere assume m 2 > 0. When m 2 < 0 
(Goldstone case) a vacuum reconstruction is needed which leads to changing both 
the initial perturbation theory and all the equations of the approximation scheme. 


i) d = 2. We will seek solutions for which the representation £(xy) = T,S(x — y) 
is valid, £ being a constant. Then equation (]34|) has the following solution in the 
momentum space: 

A = -yW+-A-r, (48) 


p 2 + p 2 m 2 + p 2 


where 


A = 


rrx 


n 


/T 


B = 


m 


m 


n 


= 1 - A 


m“ — p* 

and the notations are introduced: 


m- 




(49) 


p 2 = m 2 + £, m' 2 = m 2 + C = m 2 + E c [m). 


(50) 


The operator £ built from A in accordance with (^), (|3^) and ( |48| ) reads 

£ = AE c [p] + BE c [m] = A(E c [p] - £ c [m]) + £ c [m], (51) 


Taking into account formula ([IT]) and dehnition (|50|) leads to an equation for p 2 


A m 12 — p 2 
2 n p 2 — m 2 


In 


9 

m 



/2 



(52) 


Equation (^) has two solutions. One of them is p 2 = m' 2 which implies A = 0, 
B = 1 hence A = A c . This is the trivial solution. Another one is a solution of the 
equation: 


A 1 m 2 

27r p 2 — m 2 p 2 


1. 


(53) 
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A solution of (|53|) exists if A is negative and at small A can be written as 

2 2 ( 27T??7 2 \ 

li ~ m exp I-I . 

The corresponding values of A and B are 


( 54 ) 


in 


/2 


A ~ —- = 1 + 0(A), B ~ 


m 


and propagator is: 


A ~ 


1 


m — m 

- 2 = 0 ( A ) 

(55) 



+ 0( A). 

(56) 


/i 2 + p 2 

The following circumstances connected with solution (|54| - ^6|) should be mentioned. 


First, the non-perturbative character of the solution is obvious. Second, the negative 
value of the A parameter does not obligatory mean that the value of the physical cou¬ 
pling constant is negative. To answer the question about the vacuum stability and 
the sign of physical coupling corresponding to this solution one should at least con¬ 
sider the next step of the approximation scheme when [f? 2 ] con 7^ 0- Such a situation 
is not unusual one in the quantum field theory (see 0, for example). 

ii) d = 3. Let us try to find a solution X(xy) = Eh(x — y) again. All formulae 
©-© remain to be valid in this case too. The only change consists in using 
formula ([43]) for the two mass operators difference. The equation for y 2 has the 
form: 

A m' 2 — y 2 


J2 


m — y = 


2ti m + y 


(57) 


Besides the trivial solution y 2 = in' 2 (A = A c ) the non-perturbative one exists: 


2n' 


y = —m 

Since y = \fy? this solution exists at A < —27 t m. 

iii) d — 4. Now £ c = —A p 2 + C where we have denoted 

A = 


(58) 


A 


16n 2 

A solution should be sought in the form X = ap 2 + b. It still looks as 

m' 2 — (1 — A) y 2 


but now 


A = 


where 


(1 + a) (in 2 — y 2 ) ’ 

2 m 2 + b 

V = -j—— , 
1 “I ~ CL 


(1 — A) m 2 — m' 2 
(1 + a) ( m 2 — y 2 ) ’ 

(59) 

m 12 = m 2 + C. 

(60) 
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It follows from definitions (|35| - p6| ) and equations ( fl8|) and (|59| ) that 

E = ap 2 + b = AYj c [p] + BY, c [m] = 

= A (^c[lA ~ ScH) + y-^ S c[™]- 
It leads to the system of equations for a and b: 

a (1 + a) — A (A — 1) 


b = AC' + \—- C. 

1 ~h CL 


( 61 ) 


(62) 


We hnd the two solutions for a: aq = —A and 02 = A — 1 . For the first of them 
A + B = 1 and p 2 = m 2 + 0(A), a = 1 + 0(A). This solution is perturbatively close 
to the free one: 

A = A c + 0(A). (63) 

The second solution gives 

A + B= ^ T X \ A=l + 0( 1). 

A A 

This solution is of the non-perturbative character: 


A 


1 

A (/i 2 + p 2 ) 


+ 0 ( 1 ). 


(64) 


Besides above mentioned solutions there can exist the ones with a dipole term: 


A 


A' 

m 2 + p 2 


B' 

(m 2 + p 2 ) 2 


(65) 


We will not dwell upon this case which as a matter of fact is the particular case of 
the previous consideration when b —> m 2 a. We only point out that the dipole term 
is always of the higher order of the smallness in comparison with the simple pole 
contribution: B' = o\(A'). 

To the end of this section we would like to discuss the regularization dependence 
of our results. First, such a dependence is obviously absent when d < 4. Indeed, 
introducing quantities 

A' = A - A c , E' = E — E c , (66) 


in place of A and E, we obtain from (p6|) the following equation: 


A' = -(m 2 — d 2 + E' + Ec)" 1 * E' * A c . 


(67) 


One can see from ( |67| ) that at d < 4 the quantity A'(x = 0) is a finite one and 
therefore so is E'. Having solved equation O we get, of course, all the same results 
as we do above, irrespective of the regularization. 


14 











In case of d = 4 the situation is more complicated. When using the standard 
regularization schemes like cutoff one we find qualitatively the same results but the 
problem of cutoff removing arises. In our approach this problem do not appear at 

all. 

Conclusion 

The main idea of the proposed approach is to enlarge the class of solutions of 
the Schwinger equation by widening the set of admissible boundary conditions. But 
the problem of boundary conditions itself remains unsolved, of course. In essence, 
when constructing the approximation scheme in section 1 we constrained ourselves 
to the solutions which “copy” the connected structure of the perturbation theory. 
But even under such hard restriction it is possible to find the solutions of the non- 
perturbative character (section 2). An appropriate choice of the solution can be 
made in comparing it with the physical phenomena are to be described and on the 
base of general physical principles such as the minimal energy value of the ground 
state. 

In this connection the investigation of the next step of the scheme containing 
the two-particle amplitude (equations ([26]) - (|27|) ) is of undoubted interest. Also it 
would be interesting to apply this method to other more sophisticated quantum field 
models. 

At last it is worth mentioning that one can construct an approximation scheme 
taking as an input the loop expansion or the kinetic expansion |T^[ instead of the 
perturbation theory. 

The summation of the perturbative series, independently of how far ahead this 
problem could be elaborated cannot give us all the variety of solutions of the quantum 
field equations. These equations contain much more information than the sum of the 
perturbation theory does. We hope that the development of approximation schemes 
taking into account the singular character of the interacting fields equations allows 
one to enlarge the number of the phenomena described by the quantum field theory. 
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